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Abstract 

We discuss medium modifications of the unitarized pion-pion inter- 
action in the nuclear medium. We incorporate both the effects of chi- 
ral symmetry restoration and the influence of collective nuclear pionic 
modes originating from the p-wave coupling of the pion to delta-hole 
configurations. We show in particular that the dropping of the sigma 
meson mass significantly enhances the low energy structure created by 
the in-medium collective pionic modes. 

1 Introduction 

Modifications of hadrons properties in nuclear and hot hadronic matter is 
one of the central subjects of present day nuclear physics. For instance, it 
has been suggested in Q that the in-medium pion-pion interaction might 
be significantly reshaped at density even below normal nuclear matter den- 
sity. The basic mechanism is linked to the nuclear collective pionic modes 
sometimes called pisobars originating from the p-wave coupling of the pion 
to delta-holes states. According to detailed calculation |2|, ^ these collective 
modes are needed to explain charge exchange data ^, despite the periph- 
eral character of these experiments. This medium effect yields a softening 
of the pion dispersion relation and consequently a modification of the two- 
pion propagator involved in the unitarized T matrix describing pion-pion 
interaction at finite density. Indeed, on the basis of purely phenomenolog- 
ical models Q, an important reshaping of the pion-pion interaction in the 
scalar- isoscalar channel (sigma channel) , producing a sizeable accumulation 
of strength near two-pion threshold, has been predicted. This problem has 
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been reinvestigated with chiral symmetric models such as hnear or non hn- 
ear sigma models Q with special emphasis on the consistency between chiral 
symmetry constraints and unitarization. It has been soon realized that this 
medium effect is of considerable importance for the still open problem of nu- 
clear saturation since an important part of the nucleon-nucleon interaction 
comes from correlated two-pion exchange and several papers have brought 
extremely interesting results ^ |l^ . A possible evidence for this reshap- 
ing of the vr — TT strength function is provided by the vr — 27r data obtained 
on various nuclei by the CHAOS collaboration at TRIUMF [^]. Recent 
calculations show that the observed marked structure in the vr+vr" invariant 
mass spectrum can be partially explained by this reshaping |1T2| , p^ ]. These 
last results have been questioned in a recent paper |14| where it is found that 
pion absorption forces the reaction to occur at lower peripheral density. . 

However, what was ignored in the previous approaches was the possible 
medium modification of the basic vr — vr interaction i. e. the vr — vr poten- 
tial from the underlying quark substructure or, in other words, from the 
in-medium modification of hadronic properties associated to chiral symme- 
try restoration. This may have considerable consequences since, applying 
by hand in the linear sigma model pion-pion potential a Brown-Rho scaling 
of the sigma mass yields a significant enhancement of the near threshold 
structure in the 2tt strength function ||l5[. To go further, it is important to 
construct an in-medium vr — vr pion-pion potential This question has already 
been addressed [16| in the framework of the Nambu-Jona-Lasinio (NJL) 
model at finite temperature. In this paper we will first use exactly the 
same scheme but with a straightforward generalization at finite constituent 
quark density to be identified later with one third of the baryonic (i.e. nu- 
cleonic) density. In particular we will show that, using a slightly different 
prescription for the loop integrals, the NJL model regenerates the linear 
sigma model pion-pion Born term amplitude. This result, which is numeri- 
cally extremely close to the more involved calculational scheme of ref . jl6U , 
can be generalized at finite density and/or temperature provided the values 
of the pion mass, the sigma mass and the pion decay constant are replaced 
by their in-medium values calculated in the NJL model. Using this scheme 
we are in position to study in a particularly simple way direct observable 
consequences of both partial chiral symmetry restoration such as dropping 
of the sigma meson mass and collective p-wave pionic modes by looking at 
the in-medium pion-pion interaction in the scalar-isoscalar channel. This is 
of utmost importance since, as emphasized recently by Hatsuda et al, the 
evolution of collective scalar-isoscalar modes, i.e. the sigma meson, may 
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reveal precursor effects associated to chiral symmetry restoration |17|. 



Prom this density dependent effective linear sigma model potential im- 
plemented with a phenomenological form factor, it is possible to construct 
a unitarized vr — vr scalar-isoscalar amplitude which both preserves chiral 
symmetry constraints (Weinberg scattering length in the chiral limit) and 
reproduces experimental phase shifts. On top of precursor effects of chiral 
symmetry restoration, the inclusion of medium effects associated with the 
modification of the pion dispersion relation should be consistently done in 
the framework of the NJL model by direct coupling to constituent quarks 
populating the Fermi sea. However, the resulting pion p-wave polarizability 
calculated with quark-particles-quark-holes would completely miss the phe- 
nomenologically well established strong screening effects from short-range 
correlations {g' parameter). Incorporation of correlation effects in the NJL 
model obviously requires a much more involved level of sophistication, hence 
loosing the simplicity which is one of its main interest. Furthermore we will 
calculate the in-medium two-pion propagator from standard pion-nucleus 
phenomenology. In other words, the p-wave pion polarizability will be taken 
as its nuclear matter expression dominated in the region of interest by the 
A-hole piece corrected by screening effects. The underlying philosophy can 
be summarized in saying that the medium modified soft physics linked to 
chiral symmetry (mj^ , fj^ , low energy tt — vr potential ) is calculated within 
the NJL model while p-wave physics yielding pionic nuclear collective modes 
is described through standard nuclear phenomenology. 

2 The TT — TT Interaction from the NJL Model 

To examine the vr — vr interaction at finite density (and finite temperature). 



we start from the SU(2) version of the well-known NJL model|18]: 

C = i;{i^- mo)V' + oKi'i^f + (ipiTorijf] , (1) 

where 5 is a coupling strength of dimension [Mass]~^, and mo is the cur- 
rent quark mass. To maintain the asymptotic freedom at high energies, we 
simulate the situation by regarding the coupling constant in Eq. (|l|) as a 
momentum dependent one g(p): 

4 

9{p)=9Y[(^{A-M , (2) 

i=l 
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where Pi's are the momenta of quarks and 9 is the step function. Three 
parameters, namely g, A, and mo are determined to reproduce the pion 
mass rriTr, the pion decay constant 7,^, and the mass of sigma meson m^j. 
Here for and fj^ we choose the well- accepted values, and for we 
select two different empirical values to test the influence of the input data 
on our final results. In Tab. ||, we list two sets of values for m,r) /tt and mo-, 
and also (7, A and mo which are fixed respectively to give these two sets of 
values of m^^, f-w and m^j. In the parameter fitting, we follow the formalism 



in the paper of Hatsuda et al. |1S], where the three-momentum cut-off was 



used. With the lowest value of nia = 700 MeF, one has a set of parameters 



very close to the one of Quack et al. |16| used for finite temperature studies. 



In the following mainly devoted to finite density, we will prefer the second 
set with rria = 1 GeV for two reasons. On one hand, this second set gives a 
higher density for chiral symmetry restoration (about 2.6/9o in place of 1.5/?o) 
and on the other hand, the values of mo- = 1 GeV has been successfully used 
in previous works to fit the vacuum vr — vr data 0, |l|] 



In Tab. |l|, we also include the results with different schemes of calcula- 
tion, the "exact" one and the "Simplified" one. For the "exact" scheme, we 



just do the same calculation as Quack et al. [16| but at finite density, while 
for the "Simplified" one, we neglect the g^-dependence in all of the integrals 
(for details see Sec. 3) in our calculation, where q is the external momentum 
of pion and sigma mesons. 

The main purpose of this section is to get the density-dependence on the 
vr-vr scattering length, a^, which is related to the vr-vr scattering amplitude 

at threshold as: 

= i — ReT^ (3) 

327rm,r 

where I is the total isospin. 

To the lowest order in l/N^ {Nc, the number of quark colors), the in- 
variant amplitude Tabbed of the vr — vr scattering process (a, b and c, d are 
the isospin indices) is calculated from the box and cr-propagation diagrams 
shown in Fig. |^ (see Fig. 1 in Ref. [|l^). Following the notation of Ref. 



[iBI, we have 



%b;cd =< cpc, dpd\T\apa; bpb >= A{s, t, u)SabScd+B{s, t, u)SacSbd+C{s, t, u)5adhc 

(4) 

where s, t and u are the usual Mandelstam variables: s = {pa -^-Pb)"^-, t = 
{Va - Pc)^ and u = {pa- Pdf ■ 

Due to the crossing relations, there are three possibilities for the box 
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diaffram pO] which, aft er a direct evaluation, gives: 

{Ti)ab;cd = -{^ab^cd + &achd " ^adhc) [AN.Nfigl^^] [/(O) + I{p) - p^K{p)] 

{T2)ab;cd = -{^abScd " ^achd + 5adhc)[^N,N figlq^][I{{)) + I{p) - p^K{p)] (5) 

{r3)ab;cd = -i-SabScd + Sac^bd + 5ad<5fe,) [8iV,iV/iff^,,] [1(0) +p^L{p)/2 - 2p^K{p)] 

in terms of integrals /(O), /(p), K{p) and L{p) (see appendix). 

The diagram with intermediate a propagation shown in Fig. |^ can be 
simply expressed in terms of its components. By combination, one has : 

{%)ab;cd = -8abScdgUcS^''''^iP^-P)fD^{2p) 

{%)ab,cd = -{5achd + Kdhc)gUq[T''''^{p,p)?D„{Q) (6) 

where 

V^^ip^-p) = -8NcNfmI{p) 
r---(p,p) = -8N,Nfm[I{0)-p^K{p)] (7) 

are the (T-tt-tt vertex in the s-channel and ^-channel respectively, and 

■J 

^''^^^ " 2N^Nf[{k^ - 4m2)/(A:) - mll{m^)] ' 

is the sigma meson propagator. Similar to Eq. (8), we also have ]l6| , |20[ ] 

g;^\ = -N^[IiO)+Iip)-mlKip)f, with (N = N,Nf) . (9) 

We note here that all above results are evaluated at threshold {pf = p^ = 
= s/4), and we restrict them, for T > 0, to the scattering of pions whose 
cm. system is at rest in the heat bath|16|. 

Summing Eqs. (5) and (6)and with the help of projection of the ampli- 
tudes on total isospin /, one finds the s-wave scattering amplitudes: 

T° = 6T3 + 3T4 + 2T5 
= 

T2 = 2T3 + 2T5 (10) 

where the % are the function in Eqs. (5-6) excepted of the isospin factors. To 
get the density dependence on 7^, we must calculate the constituent quark 
mass m and the pion mass m,r at finite density first. 
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As in Ref. |1£], we get the in-medium constituent quark mass from the 
gap-equation due to the one-loop quark self-energy diagram (Fig. 5 in ref. 
[[19|). The in-medium pion mass (and also cj-meson mass) is determined 
from the dispersion relation (in the g — > limit) of the meson excitation 
(which itself is the solution of the Dyson equation in the ring-diagram ap- 
proximation). For the details of deriving the gap-equation for quarks and 
the dispersion relation for vr (and also cr) meson, one can resort to Ref. 



[16, ^ with a straightforward extension at finite density. Some specific 
results and integrals are explicitely given in the appendix . Here we display 
only the final results evolving with respect to density at zero temperature. 
One thing to be mentioned here is that there are two different definitions for 
the coupling between pions and quarks, namely g-wqq'- one is from Hatsuda 
et al. 1 19 1 and the other is from Quack et al. [|^. These two g-j^qq differ 



from each other only of a term proportional to 

/? / dk^f{E)f{-E)F{E,p) 

where f3 is the inverse of the temperature, f{E) the Fermi-distribution func- 
tion and F{E,p) a rational fraction oi E {E = \/ m? + fc^) and p (external 
momentum of pions). It turns out that this term disappears when the tem- 
perature goes to zero. So, the discrepancy between these two definitions will 
not give rise to any difference in the parameter fitting which is done only 
at zero temperature and also in the main part of this paper as we focus our 
calculation on the cases of finite density but at zero temperature. 

The density dependence of constituent quark mass, vr mass, a mass and 
/,r is shown in Fig. |2[ We find that all these quantities evolve with density 
almost linearly except in the high density region (say p > 0.30 fm^'^) and the 
trends of these curves are all the same as those evolving with temperature 
[19|. Here we use two sets of parameters listed in Tab. [l|. 

We first show that the density dependences of the threshold amplitudes 
and scattering lengths are very similar to those obtained at finite temper- 
ature in |T^. In Fig. ^, we draw the curves of the various pieces of the 
scattering amplitudes as a function of density (at zero-temperature) with 
both sets of parameters. Again, we see that all the curves evolve linearly 
with density in the low and medium density region. At a certain density pd, 



there exists a divergence as at high temperature |16|. The similar behavior 
also holds in Fig. ^ where the relation between scattering length and density 
is illustrated. An interesting feature is that one gets a higher value of density 
Pd with the higher input vacuum cr-meson mass. This divergence occurring 
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in the sigma propagator actually corresponds to the point where the in- 
medium sigma mass drops to twice the in- medium pion mass. In the second 
set of parameters the density pd turns out to be pd ~ 2ApQ. Once unitariza- 
tion is done (see section 5) this density goes down to about pd ~ 2po due to 
the presence of the two-pion loop in the sigma propagator. We will come to 
this important point in section 5 since, as pointed out by Hatsuda et al. |17], 
the vanishing of the real part of the inverse sigma propagator at ii^ = 2771,^ 
can be seen as a precursor effect of chiral symmetry restoration. However 
the aim of this paper is not primarily to study the phase transition region 
but the moderate density region below nuclear matter density accessible in 
the TT — 27r experiment. For that purpose, we show in Fig. 5 the low density 
behaviour of the scattering lengths. 



3 The Simplified Scheme 

Now we turn to the simplified scheme of our calculation, where we neglect 
the momentum dependence in all of above integrals, i.e. I{p), K(p) and L{p) 
(we will state in the later of this section the reason for such a treatment). 
But before this simplified processing, we would like to check whether the 
Weinberg limit is still kept at finite density and finite temperature. First we 
want to see how these integrals (and then other related quantities) behave in 
the first order chiral expansion (expanding to the order of m^). Keeping in 
mind that we are always in the 3-dimensional cut-off constraint and following 
the regularization of Ref. ||2^, we find after straightforward calculation, that 

I{p^=ml) = I{0)-^mlK{0) , 

Kip' = ml) = KiO) - ^mlm , 

L{p'=ml) = m-^mlh2{0) (11) 

where 132(0) is shown explicitly in the appendix. With these expansions, we 
get from Eqs. (5-9) : 



A^[2/(0) - lmlK{0) - ml{K{0) - ^mlL{0))] ' 
^2 = Ti , 

^ ^ 8i[/(0) + miim - |m^J32(0))/2 - 2m^(i^(0) - |m^L(0))] 
' iV[2I(0) - lmlK{0) - mUK{0) - im2L(0))] 
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-8m(/(0) - \mlK{f)))f 



[27(0) - \mlK{Q) - mliKiO) - '^mimW 



2' 

i 



2iV[(4m2 - 4m2)(I(0) - ^mlK{0)) - (1(0) - \mlK{0))] 



3 

-8m(/(0) - m2 (K(0) - im2L(0))]2 



[21(0) - lmlK{0) - ml{K{0) - ^mlL{0))]^ 
i 

2A^[-4m2l(0) - m2 (7(0) - ^mlK{0))] ' 



(12) 



Keeping to the first order of with the help of Eqs. (3) and (10), one 
finds finally that 



1 7ml 



327rm^ /2 

a' = (13) 



327rm^ /2 



where /tt is defined as in Ref. [19|. Note there that the function 7(0), K{<d), 
L(0) and 732(0) are the integrals at finite density and finite temperature. 
We thus recover the Weinberg limit (in the 1st order chiral expansion) at 
threshold for finite density and finite temperature. 

In the simplified calculation, we neglect the momentum dependence of 
all above integrals inside the formula for %, i.e., = = except that 
we keep k'^ in Eq. (8) unchanged as the momentum transfer in the two-pion 
interaction. Thus one gets 

^ ^ ^ 2i 
Ti = T2 = T3 



iV7(0) 
8im^ 1 

- - Nm fc2 - ml ^^^> 

where = Am? + in the simplified calculation, which can be easily 
verified from Ref. [jl^. Since /2 = —2iNm'^I{0) in the simplified calculation, 
the above amplitudes can be rewritten as : 

% = %^i\'fi ^ ^ (15) 
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with the CTTTTr couphng constant A defined by 2A/^ = — as in the 
hnear sigma model (see for instance sect. 4 of |0]). In other words, this 
simpUfied scheme yields exactly the amplitudes of the linear sigma model 
but with A, rria and m^r being functions of density (and temperature). 
The numerical results for the simplified calculation are also displayed in 
Fig. ^ and § by dashed lines. We find that these two calculations have 
no much difference from each other except for the high density region. So, 
the simplified calculation is a rather good approximation to address in our 
region of interest the density dependence of m, m^r, m^, f^^ and even the 
scattering amplitude and scattering length. 

4 In-medium Unitarized T Matrix 

We now incorporate the previous vr— vr amplitudes considered as a irreducible 
potential into a Lippman-Schwinger (LS) equation to get the unitarized T 
matrix in the scalar-isoscalar channel. However, when solving the LS integral 
equation, the use of the "exact" scheme with its full momentum dependence 
turns out to be of hopeless complexity. For this reason we prefer to use the 
potential from the "simplified scheme" which has been shown to be very 
close to the observables (pion decay constant, masses, scattering lengths) 
and their density dependences calculated with the "exact" scheme. Hence, 
the potential will be simply the linear sigma potential. For simplicity we 
restrict ourselves to an interacting pion pair with zero total momentum. In 
the scalar-isoscalar channel, the potential reads with obvious notations : 

< k, -k|V(E)|k', -k' >= .ik).(k')'^^^ { 3^ + + ^) 

(16) 

The prescription for the Mandelstam variables to be used in the LS equation 
are the ones used in section 4 of . E = y/s is the energy variable at which 
the unitarized T matrix is calculated, t and u are chosen according to their 
on-shell values : 

t = 2ml - '^^k^k' + 2k • k', u = 2ml - "^^k^^k' - 2k • k' (17) 

The form factor v{k), which is assumed to simulate the momentum de- 
pendence of the vertices, is to be fitted on experimental phase shifts and 
scattering lengths in the free case. The unitarized T matrix is the solution 



9 



of the LS equation 



< k, -k|r(^)|k', -k' >=< k, -k|i^(^)|k', -k' > 

+ \J < k, -k|y|q, -q > G2AE, q) < q, -q|^|k', -k' >(18) 

where the two-pion propagator in vacuum is [Q] : 

G2AE, q) = - p2 A 2^- (19) 

Ignoring the t and u dependence in the denominators, the potential (Eq. 
(|l6|)) reduces to a combination of separable potentials. Taking parame- 
ter set II (rrio- = 1 GeV) implemented with a one-parameter form factor 
v{k) = 1/(1 + k'^/GAm^), one exactly recovers the potential used in [^. 
The phase shifts are correctly reproduced up to 800 MeV and the scatter- 
ing length, ao = 0.23m~^, is also in agreement with experimental data. It 
can be checked that, with such a separable potential, the correction to the 
scattering length from the unitarization procedure, is of higher order in m^, 
thus preserving the chiral symmetry result (Weinberg scattering length) in 
the chiral limit. 

Prom the underlying NJL quark model, the above results for the T ma- 
trix, Eq. (|^), (1?), ([T8|), (|l9|) can be generalized at finite density by simply 



replacing the pion mass, the sigma mass and the pion decay constant by 
their in-medium values. With parameter set II, the pion mass does not 
change very much. Starting from the vacuum value niT^ = 139 MeV, one 
finds at half nuclear matter density m-,^{pQ/'2) = 140.3 MeV, and at satura- 
tion density m-,^[pQ) = 143 MeV. The variation of the pion mass is usually 
expressed in term of the s-wave pion optical potential or in term of an ef- 
fective scattering length Am^ = 2mT^Vopt = — 47r {bo)eff P- Here we find 
(6o)gjj ~ —0.01 Tn~^ which coincides with the experimental pion-nucleon 
scattering length. In reality the pion optical potential is more repulsive due 
to higher order rescattering effects but still giving very small modification to 
the pion mass. Hence, this modification of the pion mass, which simulates 
s-wave pion-nucleus interaction is practically negligible. In other words, the 
model does not give unrealistically large shift in the pion mass. 

On the contrary, there is a sizeable dropping of the sigma mass of 110 
MeV at half nuclear matter density and 225 MeV are normal density. Con- 
sequently, one finds a significant reshaping of the vr — vr strength function, as 
we will see in the next section, even at p = O.bpo which is the typical density 
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reached in vr — 27r experiment [13|. We now come to study the dressing of 



the pion propagator by its p-wave coupUng in matter. As already stated in 
the introduction, the direct use of the standard NJL model for this problem 
would give totally unrealistic results since there is no repulsion and short 
range correlation {g' parameter). Hence, we turn to a conventional nuclear 
matter description of the p-wave self-interaction of the pion on top of the 
NJL approach for the basic pion-pion potential. The pion propagator has 
the form : 

D^{\l,uj) = [uj"^ - ml -]^^ - S^{\l,uj)\-^ (20) 

where S'7r(k, w) is the conventional nuclear matter p-wave pion self-energy to 
be specified below. Notice that m.,^ is the in-medium pion mass calculated in 
the NJL model. According to the previous discussion, the small s-wave pion 
self-energy is partially accounted for but this gives in practice a negligible 
effect for the in-medium vr— vr T matrix. For this first work, mixing the effects 
from the underlying quark structure (modification of ma,TT,fTr) with those 
from pion p-wave coupling, we use the simplest possible approach. Hence, 
in the spirit of Ref. Q , we assume that the pion self-energy is dominated 
by virtual A — h excitations in the domain of energy of interest, namely : 

s^ik, uj) = e oj) = e nO(k, u:)/ ii- 5AAn°(k, u)) (21) 



Here ~ 0.5 accounts for the short range screening of the A — h polariza- 
tion bubble H*'. The numerical calculations are performed in the framework 
of the two- level model already used in Ref. The main approximation 
of this model is to neglect the Fermi motion of the nucleons in the A — h 
bubble. In addition we first neglect the width of the delta resonance. This 
last approximation will be relaxed below. The polarization bubble has the 
form : 

tf(k,^) = 1 ( p ( — i— ) (22) 



with e^k = Y + ~ ^-^^ ^{^) is the ttNN form factor. It is a 
simple matter to show that the pion propagator takes the very simple form 
typical of a two-level model : 

DAK p) = , + . %fr''\^' 

— w; p) + IT] uj'^ — (k, p) + irj 

where the eigenenergies ili, and the strength factors Zi, Z2 {Z\ + Z2 = 1) 
are very simple functions of the density which can be found in Ref. H. As 
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explained in details in this last reference, the pion strength function splits 
into two collective eigenmodes. Of special interest is the lower branch (fii), 
sometimes called the pionic branch responsible for the softening of the pion 
dispersion curve. To calculate the medium modified vr — vr T matrix, we have 
to solve the LS equation (Eq. (18)) with a modified two-pion propagator 
according to : 

^) = / ^ ^0) ^.(-k, E - ko) (24) 

In the two-level model, its explicit expression is : 

A n,{k) + n,{k) z,{k)z,{k) 



Obviously, in the absence of p-wave coupling, one recovers Eq. (19). In a 
more realistic description the main features of the two-level model survive 
but the modes fli^2 acquire a width. This can be incorporated through the 
replacement in the pion propagator (Eq. (23)) : 

+ (20) 

where Imn*^, calculated along the line j, takes into account the delta width, 
corrected from Pauli blocking, together with extra 2p-2h contributions which 
are not reducible to a delta width piece. Details are given in Ref. [21|. The 



imaginary part of the two-pion propagator is obtained through a spectral 
representation : 

1 

Im G2n (k, ^) = / du ImD^ (k, E) ImD^ {k,E-u;) (27) 

And the real part is calculated with a dispersion relation. 



5 Results 

To clearly illustrate the basic physical mechanisms, we present the results 
for the in medium pion-pion scattering matrix in the simplified calcula- 
tional scheme discussed in the previous sections. First, the treatment of 
collective vrA^A configurations entering the two-pion propagator is done in 
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the extended two-level model summarized just above. In addition we use 
the leading order term of the so-called \/N expansion of the linear sigma 
model. In a forthcoming paper we will incorporate the full linear sigma 
model potential together with a full calculation of the in-medium two-pion 
propagator. However, we believe that the conclusions concerning the rela- 
tive weight of collective pionic modes and chiral symmetry restoration will 
remain valid. 

The linear sigma model can be seen as a 0{N + 1) model with = 3. It 
has been shown |l2| that, to leading order in a expansion, one gets 
a consistent symmetry conserving approach fulfilling Ward identities and all 
chiral symmetry constraints. The corresponding potential is obtained from 
Eq. (|l6|) by simply keeping the s channel pole term and dropping the u and 
t sigma propagators. In practice it has the great advantage of making the 
potential separable. According to |T^, we take the potential as : 



<k, -k|y(^)|k', -k'>=6^;(A:)t;(A;')A5— ^ with \='\-^ 

(28) 

with a phenomenological form factor v{k) = g/{l + fc^/g^)". It is important 
to note that the introduction of such a form factor does not destroy in any 
way the chiral symmetry properties. Taking g = 0.9, a = 3 and = 1 GeV 
one gets a reasonable fit to experimental phase shifts With such a 

separable potential one easily obtains the explicit form of the T matrix : 

<k, -k|T(^)|k', -k' >=6v{k)v{k')X- E^-ml 



E'^ -ml-3X{E'^ - ml) J:{E) 

(29) 

with S(£^) given by : 

nE)=J ^v\q)G2M,E) (30) 

On figure 6, we present the results of the calculation for the vr — vr strength 
function (i.e. the imaginary part of the scattering amplitude) for zero den- 
sity, half nuclear matter density (0.5po)a'nd normal nuclear matter density 
{po)- We limit ourselves to the low energy sector {E < 600MeV) since it 
corresponds to the domain of validity of the low energy effective theory and 
to the region experimentally probed for instance in the CHAOS experiment. 
On Fig. 6a we keep the sigma mass, the pion mass and the pion decay 
constant at their vacuum values and incorporate the effects of p-wave ttNA 
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collective modes. We recover the well-known structure near the two-pion 
threshold originating from the softening of the pion dispersion relation and 
extensively discussed in [Q, ^, 0- In Fig. 6b, we disregard the effect of these 
pionic modes but simply replace in the vacuum T matrix the sigma mass, 
the pion mass and the pion decay constant by their in-medium value. In this 
way we isolate the effects intimately related to chiral symmetry restoration. 
We see a spectacular enhancement of the strength with increasing density 
which is roughly uniform in the energy domain considered. The origin of 
this effect is mainly the dropping of the sigma mass. However, it is less pro- 
nounced in the density domain typical of ordinary nuclei than in the work 



of Hatsuda et al. |17] for at least two obvious reasons. First, we start with 
a vacuum sigma mass of 1 GeV (used in our previous works to fit the phase 
shifts) significantly larger than the values of 550 MeV or 750 MeV used in 
[p!7| . Second, the cttttt coupling constant A = m^ — slightly decreases 
at variance with the work of Hatsuda et al. where it was density indepen- 
dent. We will come to a more detailed comparison of the two works when 
discussing the sigma propagator itself. In fig. 6c, we include simultaneously 
the two mechanisms. We recover the typical low energy structure associated 
to pionic collective modes but significantly reenforced by the effect of chiral 
symmetry restoration. 

To study the evolution with density of scalar-isoscalar modes, it is also 
very interesting to study the sigma meson spectral function i.e. the imag- 
inary part of the sigma propagator. It is a simple matter to obtain the 



explicit form of the sigma meson propagator in this simple model [12| : 



with the sigma self-energy given by : 

^'^(^^ - l-3AS(i?) 

The results of the calculation are shown on figure 7 for the three previous 
cases : pionic collective modes only (Fig. 6a), chiral symmetry restoration 
only (Fig. 6b) and both (Fig. 6c). It is also interesting to compare our 
results with the calculation of Ref. where chiral symmetry restoration 
is implemented by only dropping the sigma mass, ignoring the modification 
of the fTvrvr coupling constant. Since we find in this present work that A is 
only slightly modified (about 20% lower at p = po) it is not surprising that 
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the reinforcement of the low energy enhancement is quaUtatively similar in 
these two works at least at normal nuclear matter density. 

Until now we have considered the density regime accessible in experi- 
ments on ordinary nuclei such as the CHAOS experiments. It is however 
extremely interesting to go beyond, although the basics nuclear physics in- 
gredients of the calculation are less under control. We now follow for a while 



the arguments of [17|. Let us consider the spectral function of the sigma 
meson : 

1 1 ImS„{E) 

(33) 

where rn-o- is as usual the in-medium sigma meson mass obtained here from 



the NJL model. As in the work of Hatsuda et al. |17], this in-medium sigma 
meson mass has an evolution which follows linearly at low density the scalar 
condensate. However the precise law of evolution is different since it comes 



from a nucleonic tadpole diagram in |17]. Near the two-pion threshold the 



phase space factor yields to the well-known behavior : 

ImD„ ~ ImS„ ~ {E'^ - 4m^)^/2 (34) 

Before chiral symmetry restoration with complete a — tt degeneracy, it must 
exist a density at which ReD^^{E = 2m^) = or 4m^ — — S{E = 
2m,r) = 0. At such a density pd, ignoring the effect of pionic collective 
modes, the spectral function behaves as : 

. 2,„.) = -i ~ HE - 2..,)HE- - 4™5)'/^ (35) 

This implies that there arises a mild integrable singularity just above the 
two-pion threshold in the medium. With our set of parameters this precursor 
effect of chiral symmetry restoration occurs at a density slightly larger than 



2pQ at variance with |17] where this density was found at 1.25/9o- In addition 
when collective pionic modes are included, there is strength below two-pion 
threshold due to the various sources of width (A — h, 2p — 2h). As a conse- 
quence, the mathematical singularity disappears. However, a precursor ef- 
fect still manifests itself through a strong enhancement of the near-threshold 
structure created by the p-wave pionic collective modes. This feature is il- 
lustrated on Fig. 7 where the influence of chiral symmetry restoration itself 
is visible by comparison of Fig. 7a and Fig. 7b. 
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6 Conclusion 



In this paper we have examined medium effects on the unitarized pion-pion 
interaction in the scalar-isoscalar channel and the sigma meson spectral 
function. The new feature of this work is the simultaneous treatment of 
chiral symmetry restoration and nuclear effect associated with the existence 
of collective pionic modes. Starting from the standard Nambu-Jona-Lasinio 
model generalized at finite density, we have obtained an in-medium pion- 
pion potential. In particular, we have shown that the NJL model regenerates 
the linear sigma model potential but with modified parameters (m^, and 
/tt). This modification of the parameters intimately related to chiral symme- 
try restoration strongly reenforces the vrvr strength function in the threshold 
region. Since the NJL model is not adapted to the calculation of collective 
pionic modes, we have used the standard phenomenological nuclear physics 
approach. When both effects are incorporated, the low-energy tttt structure 
resulting from the existence of nuclear collective modes is significantly en- 
hanced in particular through the dropping of the sigma mass. This feature 
is obviously of considerable importance for the understanding of the struc- 
ture observed by the CHAOS collaboration in tt — 27r experiments on various 
nuclei. 

Acknowledgments : We are indebted to Z. Aouissat, J. Delorme, M. 
Ericson, J. Marteau, P. Schuck and J. Wambach for stimulating and enlight- 
ening discussions. 
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Appendix 



We list below the various integrals appearing in the text : 



I{p) 
K{p) 

Hp)-- 

ran ip) 





1 


(27r)4 [A;2 


— m2][(A; + pY — rn?] 


d^/fc 


1 


(27r)4 [A;2 


— m?'Y[{k +p)2 — m?] 




1 


(27r)4 [A;2 - 


- rin?Y\{k +p)2 — m?Y 


d^A; 


1 



(36) 
(37) 
(38) 
(39) 



(27r)4 [/c2 - m2]™[(A; + pf - m?Y ' 

The above equations are strictly valid at T = 0. The explicit expressions 
at finite temperature are given in [||] and can be straightforwardly extended 
at both finite density and temperature. Here we simply quote a few relevant 
results for numerical calculation at finite density : 



-iK{po) 
-iL{po) 



^ dk A;2 



-1 



27r2 Evl- 4£;2 



^ dk -A;2 
p 2^4E3 
^ dk A;2 
'kp^AE^pl 



1 



8*^2 



pI - 4^2 (p2 _ 4^2)2 

1 12£;2 



(40) 

and the gap equation is : 



64£'^ 



pI - 4^2 {pI - 4^2)2 {pi - 4E^)3 



M = m + ANcNfg 



A 1,2 



k^dkM 
27r2 'W 



(41) 
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Table 1: Physical quantities at p = in NJL model for different parameters 
(I and II) and different schemes of calculations (Exact and Simplified). 

(MeV) (MeV) (MeV) (MeV) (MeV) (MeV) (fm^) 





u 






A 


m 


mo 


9 


I(Exact) 


93.0 


139.0 


700.0 


620.84 


347.49 


5.677 


0.2247 


II(Exact) 


93.0 


139.0 


1000.0 


573.31 


498.85 


5.772 


0.3303 


I (Simplified) 


93.0 


139.0 


700.0 


624.25 


343.03 


5.540 


0.2207 


Il(Simplified) 


93.0 


139.0 


1000.0 


573.60 


495.15 


5.712 


0.3284 



Figure 1: Box and a-propagation diagrams for tt-tt scattering {a,b,c and d 
are the isospin indices. 

Figure 2: Constituent quark mass m, pion mass m,r; u-meson mass mo- and 
pion decay constant vs. density (at zero-temperature) with parameter 
set I and II: solid lines for the exact calculation and dashed lines for the 
simplified calculation. 

Figure 3: Scattering amplitudes % vs. density (at zero-temperature) with 
parameter set I and II: solid lines for the exact calculation and dashed lines 
for the simplified calculation. 

Figure 4: Scattering lengths and vs. density (at zero-temperature) 
with parameter set I and II: solid lines for the exact calculation and dashed 
lines for the simplified calculation. 

Figure 5: Same as Fig. 4 but for moderate range of density. 

Figure 6: Strength function for the vr — vr interaction at zero density (dot- 
dashed curves), half nuclear matter density (dashed curves) and normal 
nuclear matter density (solid curves) . Fig. 6a (upper left panel) corresponds 
to a calculation with coupling to the pionic collective modes, keeping nia, 
m-^ and /^r at their vacuum values. Fig. 6b (upper right panel) incorporates 
the in-medium modification of m^-, m7r and f^r from the NJL model, ignoring 
the effect of collective pionic modes. Fig. 6c (lower panel) incorporates both 
effects. 



19 



Figure 7: Spectral function for the sigma meson at various densities. Fig. 
7a (upper panel) corresponds to a calculation with coupling to the pionic 
collective modes, keeping rua, m^r and /^r at their vacuum values. Fig. 7b 
(lower panel) incorporates the in- medium modification of rria, uIt^ and /tt 
from the NJL model on top of collective pionic modes. 
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